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Abstract 

We study nuclear u~ — e~ conversion in the general framework of ef- 
fective Lagrangian approach without referring to any specific realization of 
the physics beyond the standard model (SM) responsible for lepton flavor 
violation {Ifit ). All the possible types of short range interactions (non- 
photonic mechanisms), i.e. (pseudo-)scalar, (axial-)vector and tensor, are 
included in our formalism. We show that the u~ — e~ conversion in the 
strange nucleon sea via the scalar interactions is comparable with that in 
the valence quarks. This provides an insight into the strange quark cou- 
plings beyond the SM. From the available experimental data on u~ — e~ 
conversion and expected sensitivities of planned experiments we derived 
upper bounds on the generic I/it - parameters of u~ — e~ conversion sensi- 
tive to the relevant u-,d- and s-quark couplings. 
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The muon-flavor violating processes 



fj,- + (A,Z) — e" + (A,Z)*, (1) 

i.e. muon-to-electron (fi~ — e~) conversion in nuclei, is known as a very sensitive 
probe of lepton flavor violation ) and related physics beyond the standard 
model (SM) [l|]-||. This fact has been recently strengthened by the evidence 
for the muon- neutrino oscillations, drawn by the Superkamiokande experiment, 
which is the first convincing signal of the non-standard physics connected to the 
lepton flavor non-conservation. The distinct feature of coherent enhancement 
in nuclear — e~conversion makes it more promising probe of fyf than other 
lepton flavor violating processes — > e~7, etc.). In general the structure of 
a participating nucleus brings some uncertainties into the theoretical predictions 
for jjT — e~conversion. However, in the most interesting case of coherent pT — 
^conversion |7| these uncertainties can be significantly reduced by the possibility 
of using the available experimental data on nuclear densities. 

On the experimental side, at present, there is one running \i~ — e~conversion 
experiment, SINDRUM II ||, and two planned experiments, MECO ||, 10] and 



PRIME JTTJ. The SINDRUM II experiment at PSI § with 48 Ti as stopping 
target has established the best upper bound on the branching ratio 



The MECO experiment with Al is going to start soon at Brookhaven JT0|. The 
sensitivity of this experiment is expected to reach the limit 



T(u- + 27 Al -> e~ + 27 Al) .. 



This year the PSI experiment is running with the very heavy nucleus 197 Au aiming 



to improve by a factor of about 20-30 over the previous limit, R£ < 2.0 x 10 -11 , 



set on fi — e in 197 Au by the same experiment some years ago || |l2 |. Now the 
expected limit is 



r , , T(u- + 197 Au -> e" + 197 Au) ,., 



Very recently, a proposal for a new experiment at Tokyo (PRIME) was announced 
[IT]. It intends to utilize the 48 Ti as stopping target with the impressive expected 



sensitivity of < 10~ 18 ]TT 
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These experimental limits can put severe constraints on mechanisms of fi~~ — 
e~conversion. In the literature there have been studied various mechanisms be- 
yond the SM (see [0-0 and references therein) classified into two categories: 
photonic and non-photonic as shown in Fig. 1. Specific mechanisms from both 
categories significantly differ in many respects and, in particular, in nucleon and 
nuclear structure treatment. This is attributed to the fact that they operate at 
different distances and, therefore, involve different details of the nucleon struc- 
ture. Long-distance photonic mechanisms (Fig. 1(a)) are mediated by virtual 
photon exchange between the nucleus and the \i — e lepton current. They suggest 
that the fi~ — e~ conversion occurs in the lepton-flavor non-diagonal electromag- 
netic vertex which is presumably induced by non-standard model physics at loop 
level. The hadronic vertex is characterized in this case by ordinary electromag- 
netic nuclear form factors. Contributions to /i — e conversion via virtual photon 
exchange exist in all models which allow /j, — > e7 decay. On the other hand, short- 
distance non-photonic mechanisms (Fig. 1(b)) are described by the effective fyf 
4-fermion quark-lepton interactions which may appear after integrating out heavy 
intermediate states (W, Z, Higgs bosons, supersymmetric particles etc.). 

In this Letter we focus on the non-photonic mechanisms of [i~ — e~conversion 
(Fig. 1(b)). The generic effect of physics beyond the SM in /i~ — e~is described 
in our approach by an effective Lagrangian which includes all the possible 4- 
fermion quark-lepton interactions. Our special interest is concentrated on the 
scalar interactions which are sensitive to the heavy quark content of the nucleon. 
Also, the current fi~ — e~experiments at Brookhaven and PSI can efficiently probe 



the scalar component of the — e~conversion ||-[12|. We will show that the 
contribution to the — e~conversion rate which originates from the strange 
nucleon sea via the scalar interactions is comparable with that coming from the 
valence quarks of the nucleon. 

We derive a general representation of the fi~ — e~branching ratio in terms of 
generic Ifjf parameters of the effective 4-fermion quark-lepton fi~ — e~transition 
operators. Transforming these operators, first to the nucleon and then to the 
nuclear level, we pay special attention to the effects of nucleon and nuclear struc- 
ture. The nucleon structure is taken into account on the basis of the QCD picture 
of baryon masses and experimental data on certain hadronic parameters. For nu- 
clear structure calculations we apply the formalism described in Refs. 0-0- 
Our applications refer to the nuclei of current experimental interest, 27 Al, 48 Ti 
and 197 Au, with special attention to 197 Au which has not been previously studied 
in the context of /i~ — e~conversion. We start with the 4-fermion effective La- 
grangian describing the non-photonic \i — e conversion at the quark level which 
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(a) (b) 

Figure 1: (a) Photonic (long-distance) and (b) non-photonic(short-distance) con- 
tributions to the nuclear /x~ — e~ conversion. 

corresponds to the upper vertex of the diagram in Fig. 1(b). It can be written 
in a general Lorentz covariant form as 

rl ^ F r (<?) -A jBfi , (q) -A jB , (<?) ■ (K\ 

*-eff — ~7k [VABJfi J (q) + VABJ J (q) + Vt 3[LV J • V>) 

A,B,C,D;q 

where the summation involves A, B = {A, V}, C,D = {S, P} and q = {u, d, s}. 
The Ifjf parameters rfl depend on a concrete tyf model. The lepton and quark 
currents are = e-y^, j£ = e^ 5 fi, j s = e fi, j p = ery 5 fi, = ea^fi, 

J( q ) = q^q, J(g) = qi^q, J( q) = gq, J p q) = 975 q, Jj$ = q^ u q- 

The next step is the reformulation of the quark level Lagrangian ([5]) in terms of 
the nucleon effective fields. First we write down the nucleon level Lagrangian in a 
general Lorentz covariant form with the isospin structure of the \f~ — e~transition 
operator 



'"eff ~ c L [Jfj, \ a AB J (0) + a AB J (3) ) + J \ a CD J (0) + a CD J (3))^ 



V2 



A/ AO) tB^ , (3) T B^ , -Of AO) T D , (3) 



A,B,C,D 
(0) jpv , A3) 



+ ^(<^o) + 4%))J- (6) 

The isoscalar J( ) and isovector J( 3 ) nucleon currents are defined as J^f = 

N^r k N, J$ = N-y^nN, J^ k) = Nr k N } J( k) = N lb r k N, Jfc = Na^r k N, 

where k = 0, 3 and tq = I. In Eq. (||) we neglected derivative terms. In the 
matrix elements of yT — e~ conversion they produce a small contribution propor- 
tional to q/m p < m^/nip ~ 0.1 where q = |q| is the momentum transfer to the 
nucleon while and m p are the muon and the proton masses respectively. 
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Now we relate the coefficients a in Eq. (j^) with the "fundamental" ^/parameters 
7] of the quark level Lagrangian fl5|). Towards this end we apply the on-mass-shell 
matching condition |]T3[ 

^ F \Cl ff \^j) « (^CZf^j), (7) 

where and are the initial and final nucleon states. 

In order to solve this equation we use various relations for the matrix elements 
of the quark operators between the nucleon states 



(N\q T K q\N) = G% N) $ N T K 



Ni 



with q = {u, d, s}, N = {p, n} and K = {V, A, S, P, T}, T K = {7^, 7^75, 1, 75, a^}. 
Since the maximum momentum transfer in /i — e conversion is much smaller than 
the typical scale of nucleon structure we can safely neglect the q 2 -dependence 
of the nucleon form factors G^' N ^ and drop the weak magnetism as well as the 
induced pseudoscalar terms proportional to the small momentum transfer. 
Isospin symmetry requires that 

n {u,n) _ n (d,p) _ n d ri{d,n) _ n (u,p) _ n u ^(s,n) _ s~i(s,p) _ sys / \ 
U K — = ^Ki — = "JO U K ~ ^ K = ^JO V ) ) 

Note that the axial, pseudoscalar and tensor nucleon currents couple to the nu- 
clear spin leading, therefore, to incoherent contributions. Thus, only the vector 
and scalar nucleon form factors are needed in the case of the coherent nuclear 
pT — e~ conversion which we are studying in the present paper. 

Conservation of vector current implies that the vector charge is equal to the 
number of the valence quarks of the nucleon and, therefore, 

G v = 2, G d v = l, G s v = 0. (10) 

The scalar form factors can be extracted from the baryon octet B mass 
spectrum Mb in combination with the data on the pion-nucleon sigma term 
a n N = (1/2) (m u + md)(p\uu + dd\p). We follow the QCD picture of the baryon 
masses based on the relation (B\8£\B) = MbBB and on the well known repre- 
sentation [|14]] for the trace of the energy-momentum tensor 9£ = m u uu + m^dd + 

m s ss — (ba, s /d> , K)G a G% u , where G a is the gluon field strength, a s is the QCD 
coupling constant and b is the reduced Gell-Mann-Low function with the heavy 
quark contribution subtracted. Using these relations in combination with SU(3) 



relations |15fl for the matrix elements (B\9^\B) as well as the experimental data 



on M B and u^n we derive 

(7^5.1, G|w4.3, G%^2.h (11) 
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for the conventional values of the current quark masses: m u = 4.2 MeV, = 7.5 
MeV, m s = 150 MeV. These masses, however, are not yet well determined | 16| 



and uncertainty factor < 2 may affect the results in Eq. (|TT|). For our further 
numerical estimates these uncertainties are not critical and we take the central 
values Eq. (|ll|) for the form factors G^' d,s . The approach of Refs. 1 171 ) relying on 



the results of the lattice simulations, gives for these form factors similar values 
with the same level of uncertainties. What remains true is that the strange quarks 
of the nucleon sea significantly contribute to the scalar nucleon form factor Gs- 
This result dramatically differs from the naive quark model and the MIT bag 
model where G S AS P = 0. 

Now the above Eq. (0) can be solved and the coefficients a of the nucleon 
level Lagrangian (^) can be expressed in terms of the generic J^f parameters 77 of 
the quark level effective Lagrangian Eq. (|5]) as 

4? = l(v^+4 d i)(G u v + G d v ), «^ = i(^-^)(G^-G d y ), (12) 

4°J = ^+vfkGs + G d s )+ V ? s G^ af s = ^- v ^)(G" s -G d s ), 

where I = V,A and J — S,P. 

From the Lagrangian @, following the standard procedure, we can derive the 
formula for the total fi — e conversion branching ratio. In this paper we confine 
ourselves to the coherent process which is the dominant channel of fi—e conversion 
exhausting, for the majority of experimentally interesting nuclei, more than 90% 
of the total fi~ — e~branching ratio JF|. To leading order of the non-relativistic 
reduction the coherent [i — e conversion branching ratio takes the form 



G\ p e E e (M p + M r ^ 2 
47r r(/i~ — > capture) 



where p e , E e are the outgoing electron 3-momentum and energy. We define the 
quantity 

Q = Wfy + a% 0| 2 + |a<g, + «$, 0| 2 + |4°1 + 41 0I 2 + Wfs + 41 <P? 
+2 Re{(aP v + a® 0)(4°] + af s 0)* + (a% + a% <f>)(af s + a% 0)*} . (14) 

The nuclear transition matrix elements M. Pyn entering Eq. flT3|), for the case of a 
ground state to ground state \T~ — e~transition, are defined as 



M p>n = Air / j (Per)^t,(r)p p ^(r)r 2 dr , (15) 
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(jo is the zero-order spherical Bessel function) where p p , n {i r ) are the spherically 
symmetric proton (p) and neutron (n) nuclear densities normalized to the atomic 
number Z and neutron number N, respectively, of the participating nucleus. 
$ M (r) is the space dependent part of the muon wave function. 

The quantity Q in fll4|) depends on the nuclear-target parameters through 
the factor = (A4 P — M. n )/(M. P + M. n ) . However for all the experimentally 
interesting nuclei this parameter is small. Therefore nuclear dependence of Q 
can always be neglected except very special narrow domain in the Jfjf parameter 
space where < cr 3 )0. Another important issue of the smallness of the ratio 
in Eq. flTjj) is the dominance of the isoscalar contribution associated with 
the coefficients a^°\ Given that the relative significance of the strange quark 
component, which enters only the isoscalar couplings, is enhanced. 

The nuclear matrix elements Ai P7 n, defined in Eq. fll5|) , are numerically cal- 
culated using proton densities p p from Ref. |T8j and neutron densities p n from 



Ref. [|19j whenever possible. The muon wave function $ M ( r ) was obtained by 
solving the Schrondinger equation with the Coulomb potential produced by the 
densities p p .„, taking into account the vacuum polarization corrections ||. In this 
way the matrix elements M. p , n for the nuclear targets 27 A1 and 48 Ti have been 
calculated in Ref. |J. Here we apply this approach to 197 Au. The results for 
Ai Pl n corresponding to the nuclei Al, Ti and Au are given in Table 1 where we 
also show the muon binding energy e& and the experimental total rates (1^) of 



the ordinary muon capture reaction [20 



By inserting in Eq. fll3|) the values of the nuclear matrix elements Ai p>n and 
other quantities from Table 1 and Eqs. (0)-(f|) we can derive the upper limits on 
the parameters of the effective Lagrangians fl5|), @. These limits correspond to 
the existing || or expected [pi], [11] experimental limits on the branching ratio 



R^ e . The most straightforward limit can be set on the quantity Q of Eq. (|13|) but 
its physical meaning is rather obscure. In order to obtain physically interesting 
upper bounds on the ^parameters a^°\ 7}™ we adopt the usual assumption that 
different terms in expressions ( |T2"D and ( fH| ) do not substantially compensate each 
other or, equivalently, that only one term dominates at a time. In this way we 
derived the upper limits given in Table 2. The scaling factors Ba in Table 2 are 
defined as 

(j^exp \ 1/2 / j^exp \ I/ 2 

e.i-10-13 J ' Bm= ( 2 .o To-") ' Bau = 

Multiplying corresponding column in Table 2 by Ba one can reconstruct upper 
limits on the listed parameters for the case when experimental upper bounds on 
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the branching ratio R e £P differ from those we used in our analysis. 

The constraints for rjjg in Table 2 originate from the contribution of the 
strange nucleon sea. As can be seen, they are comparable to the other pT — 
e~limits derived from the valence quarks contributions. 

The limits in Table 2 represent a general outcome of the fi~ — e~ conversion 
experiments for the Tfjj physics. These limits can be easily translated into limits 
on the parameters of specific Jfjj model predicting the \x~ — e~ conversion. This 
is achieved by adjusting the quark level effective Lagrangian of the model to the 
form of Eq. (|5]) and by identifying the effective parameters r$B with expressions 
in terms of model parameters. Then the upper bounds on T]J B from Table 2 can 
be translated to certain constraints on the model parameters present in these 
expressions. 

In conclusion, we constructed a general effective Lagrangian describing the 
non-photonic fi~ — e~ conversion and specified all the fyf parameters charac- 
terizing this process without referring to any specific model beyond the SM. It 
includes (axial-)vector, (pseudo-) scalar and tensor interactions. We derived gen- 
eral formula for the coherent \x~ — e~ conversion branching ratio in terms of the 
Ifjf parameters of the quark level effective Lagrangian. We calculated previously 
unknown nuclear matrix elements of 197 Au (the current SINDRUM II target). 
We found that the yT — e~ conversion branching ratio in the strange quark sea of 
the nucleon is comparable with that in the valence quarks. We have inferred the 
generic /i~ — e~ constraints on the ^/parameters from the existing and expected 
experimental bounds on the pT — e~ conversion rate. These results are readily 
used for derivation of the constraints on the parameters of any specific Tjjj model 
without need of the nucleon and nuclear structure calculations. This provides 
useful nuclear and particle physics inputs for the expected new data from the 
PSI, MECO, PRIME and other \i~ — e~ experiments to set bounds on the muon 
I/jf violating parameters. 

This work was supported in part by Fondecyt (Chile) under grant 8000017, by 
a Catedra Presidencial (Chile) and by RFBR (Russia) under grant 00-02-17587. 
T.S.K. would like to express his appreciation to UTFSM for hospitality. 
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Nucleus 


Pe {fm- 1 ) 


e b (MeV) 


r MC (xio 6 s- 1 ) 


M p 


M n 


27 A1 


0.531 


-0.470 


0.71 


0.047 


0.045 


48 Ti 


0.529 


-1.264 


2.60 


0.104 


0.127 


197 Au 


0.485 


-9.938 


13.07 


0.395 


0.516 



Table 1: Transition nuclear matrix elements M. v<n (in fm 3 / 2 ) of Eq. (|i~5"D and 
other useful quantities (see the text). 



Parameter 


Present limits 
48 Ti(/i - e) • B Ti 


Expected limits 
27 Al(/i - e) • B Al 


Expected limits 
197 Au(/i - e) • B Au 


Q 


2.1 x 10~ 14 


1.2 x 10~ 18 


7.8 x 10~ 15 


«(°) 

AB 


1.5 x 10~ 7 


1.1 x 10~ 9 


8.5 x 10~ 8 


hv 


1.0 x 10~ 7 


7.3 x 10- 10 


5.7 x 10~ 8 




3.0 x 10~ 8 


2.3 x lO" 10 


1.8 x 10~ 8 


Us) 

IJS 


5.8 x 10~ 8 


4.4 x 10- 10 


3.4 x 10~ 8 



Table 2: Upper bounds on the I/jj parameters (see definitions in the text) in- 
ferred from the SINDRUM II data on 48 Ti [Eq. (|2|)] as well as from the expected 
sensitivities of the current SINDRUM II [Eq. (§] and future MECO [Eq. (§)] 
experiments employing 197 Au and 27 Al as stopping targets respectively. We de- 
noted: AB = VV, AV } SS, PS; I = V, A; J = S,P. The scaling factors B { are 
defined in Eq. ([IB]). 
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